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Hidden symmetry and nonlinear paraxial atom optics. 
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A hidden symmetry of the nonlinear wave equation is exploited to analyze the propagation of 
paraxial and uniform atom-laser beams in time-independent and quadratic transverse potentials 
with cylindrical symmetry. The quality factor and the paraxial ABCD formalism are generalized 
to account exactly for mean-field interaction effects in such beams. Using an approach based on 
moments, these theoretical tools provide a very simple and yet exact picture of the interacting 
beam profile evolution. Guided atom laser experiments are discussed. This treatment addresses 
simultaneously optical and atomic beams in a unified manner, exploiting the formal analogy between 
nonlinear optics and nonlinear paraxial atom optics. 

PACS numbers: 



I. INTRODUCTION. 



The realization of quasi-continuous atom lasers, build- 
ingupon a great theoretical [l[ and experimental 0, H, 
0,13 effort, is a topic of considerable excitement: be- 
yond the challenge of producing a fully coherent atomic 
beam, such devices appear indeed as promising tools for 
various applications, ranging from precision experiments 
using atom interferometry [6| to the investigation of basic 
physical phenomena such as resonant nonlinear quan- 
tum transport [81] . While the first atom lasers involved 
a matter-wave beam outcoupled into free space [2| and 
experiencing a free fall, recent experiments have demon- 
strated the possibility to obtain a guided atomic beam 
subject to an acceleration several orders of magnitude 
smaller than gravity [1, In such beams the atomic 
flux and longitudinal velocity are well-monitored, and 
the beam linear density - i.e. the number of atoms per 
unit length - is constant on the propagation axis with 
a good approximation. Interaction effects on the prop- 
agation of the extracted atomic beam are usually con- 
sidered as negligible outside the source condensate, and 
as such are discarded from the analysis. Using a de- 
scription of the atomic beam based on moments [9| , it is 
in fact possible to account exactly for mean-field inter- 
action effects in such guided atom lasers, provided that 
the atoms propagate at a constant velocity and in suit- 
able transverse potentials. The purpose of this article is 
to expose this method, exploiting commutation relations 
reminiscent of a hidden symmetry of the 2D Schrodinger 
equation [l^, I explore the connection of previous 
results of nonlinear optics [13, [3, [13] to this symmetry 
and I extend them to nonlinear atom optics, giving a 
compact proof valid simultaneously for both optical and 
atomic beams. 

The similarity between the equations of propagation 
for light and atomic waves [Tsl . [l6| , enabling one to re- 
produce with matter waves many nonlinear optical phe- 
nomena [13], has already been used to adapt success- 
ful theoretical methods from optics to atom optics [isj . 
Concerning the linear propagation of matter waves, a 



relevant example is the ABCD matrix formalism [ig ] 
introduced by Borde for atomic clouds [13] ■ This ap- 
proach computes very efficiently the propagation of dilute 
atomic wave-packets in time-dependent quadratic poten- 
tials. Among other examples, let us mention the intro- 
duction of a paraxial wave equation and of a paraxial 
ABCD formalism for atomic beams [1, |2l|j the charac- 
terization of a multiple-mode atomic beam with a quality 
factor [U, [13, [Isj , and the treatments of atomic interac- 
tions with a lensing term [1,[23] . Although the optical ap- 
proach of guided atom laser propagation is not exclusive 
- other treatments rely on hydrodynamic equations (25j 
-, it provides a reliable and tractable characterization of 
such beams. 

In this spirit, I considered two general results of non- 
linear optics concerning the transverse profile of parax- 
ial light beams. In a uniform medium, the second-order 
moment of the transverse intensity distribution follows a 
parabolic law, even in the presence of a Kerr effect 
In a graded-index medium, the transverse width oscil- 
lates with a period independent from the strength of the 
nonlinearity [13]. These results, which attest a universal 
behavior of the beam width evolution, are indeed inti- 
mately connected to a hidden symmetry [13, [HI, [13] of 
the paraxial equation. They give a simple and yet exact 
expression of the transverse size of light beams propa- 
gating in a nonlinear medium and have had several ap- 
plications in nonlinear optics [l3, [13 • As suggested in 
Refs. dH, [13]j these properties could be also relevant for 
nonlinear atom optics. 

Here, I use them to develop a simple ABCD matrix 
formalism, suitable to address the propagation of inter- 
acting atom-laser beams satisfying the following assump- 
tions: the beam is stationary, paraxial, of constant lin- 
ear density along the propagation axis, and it propa- 
gates in a transverse potential which is time-independent, 
quadratic and of cylindrical symmetry. As a preliminary 
step, we will first analyse the free expansion of such in- 
teracting atom laser beams and discuss experimental re- 
sults [|. 
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II. EQUATIONS OF PROPAGATION FOR 
PARAXIAL LIGHT AND ATOMIC WAVES IN A 
NON-LINEAR MEDIUM. 

In this section, we show that a paraxial Ught wave in a 
transparent and isotropic nonlinear medium and an inter- 
acting paraxial atomic beam of constant velocity follow 
formally equivalent equations of propagation, which take 
the form of a 2D Schrodinger equation 



2ik 



du 



-eAT^/^i + 7|'0±| ipj_ + ek a {u)r ip± 



This similarity between the propagation of light and 
atomic waves had been pointed out, in the linear regime, 
by Borde [ll]. A formal analogy between mean- field 
atomic interactions and the Kerr effect is also discussed 
in Ref. [11. 



A. Equation of propagation for optical waves. 

We consider a paraxial and monochromatic light wave 
propagating in a non-absorptive and nonlinear medium 
of graded index (GRIN). Specifically, its linear relative 
permittivity is assumed to have a quadratic position- 
dependence of the form er{x,y,z) = e°(l — 0^{z)r'^). 
This yields a medium index ri(r, z) = no(l — 0^{z)r'^)^/'^ ^ 
quadratic in the vicinity of the propagation axis where 
r <C I3~^{z). The index is greater on the propagation 
axis so that the medium acts as a converging lens and 
traps rays around this axis. The electric field is assumed 
monochromatic, written in the form 

E(a;,y,z,<) = ^i(x,2;,z)e^('=^-'^*)a, 

and it satisfies the paraxial approximation, i.e. 
|9'0j^/(9z| ^ Using this approximation in 

Maxwell's equation, one obtains the nonlinear paraxial 
equation for optical waves 28] 



2ik 



dz 



AtiP± + 7|V'-l|Vj_ - k^P^r^^^ . (1) 



The wave-vector k satisfies the dispersion relation fc^ = 
e^uj'^/c^, and the nonlinear coefficient 7 — xi^^J^ l<? is 
proportional to the nonlinear susceptibility Xe at the 
frequency to. is the transverse Laplacian acting on the 
coordinates x, y. This equation appears as a 2D nonlinear 
Schrodinger equation in which the coordinate z along the 
propagation axis plays the role usually devoted to the 
time in quantum mechanics |15|. 



B. Equation of propagation for atomic waves. 

We give a reminder on the propagation equation for a 
paraxial and monochromatic atom laser evolving in the 
stationary regime within an external potential [2]| . 



The propagation equation is derived in the general case, 
but the results presented later apply only in the absence 
of a longitudinal potential. 

Let us consider a wave- function -0, solution of the time- 
independent Schrodinger equation 

The Hamiltonian Ji accounts for interactions treated 
in the mean-field approximation and for an external 
potential C/, sum of a transverse cylindrical potential 
C/^(r, z) = m/2 x uP\jyZ) and of a longitudinal one 
C///(z). As in Ref. [1^, the wave-function is factorized in 
a 2D-I-1D decomposition as 

i'{x,y,z) = tp±{x,y,z)tp//{z). 

In our conventions, the transverse wave- function is nor- 
malized to unity Vz Jdxdy\il^±\'^{x,y,z) = 1. The 
longitudinal wave-function verifies the ID time- 
independent Schrodinger equation 



2m dz^ 

which can be solved with the WKB method as 



iP//{z) 



p{z) 



exp 



dup{u) 



T is the matter-wave flux through any transverse plane, 
p(z) — \^ 2m {E — U i/{z)) is the classical momentum 
along z, and the integral in the exponential argument 
starts at an arbitrary coordinate zg of the propagation 
axis corresponding to a point within the considered prop- 
agation zone. 

For a beam propagating in the classically allowed re- 
gion where U//{z) < E, such WKB approach is in- 
deed valid only if the wave- vector k{z) — p{z)/h sat- 
isfies the condition \dk/dz\/k <C k, otherwise quan- 
tum reflections may occur. Equivalently, the longitudi- 
nal potential U// should vary smoothly enough to verify 

\dU///dz\ < V^iE - U//)^^^/h. 

In a paraxial beam, the average momenta {p±x,y,z) 
of the transverse wave-function ip± must satisfy 
\{p±x,y,z)\ <C piz). The second derivative of the wave- 
function ij}± on the coordinate z is then negligible, and 
the wave-function tp^ satisfies 



ih 



p{z) d 



m 



dz 



A 

— At 

2m 



m 



-IV'/zH^-Ll'-c/i ^±-0 

(2) 

with fls the s-wave atomic scattering length. 

At this point, we note a significant difference between 
the nonlinear optical and atom-optical propagations. In 
typical experimental conditions in optics, the refraction 
index and the longitudinal wave-vector of the beam are 
constant along the propagation axis with a good approx- 
imation. In contrast, atom-laser experiments can involve 
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a longitudinal potential which may change significantly 
the longitudinal atomic momentum during the propaga- 
tion. This induces a difference between the optical prop- 
agation equation ([T]) and the atomic propagation equa- 
tion ([2]): in the latter, the coefficients of the first-order 
derivative and of the nonlinear term become z-dependent 
in presence of a longitudinal potential. One can elim- 
inate this dependence in the first-order derivative term 
by performing the following variable change, replacing 
the longitudinal coordinate z by the parameter 

t{z) = / dz^^ , 

defined as the time needed classically to propagate from 
the axis point of coordinate zo to the point of coordinate 
z. Eq.(I2) takes then the form 

2ik— — — ~At'4'± + 8'KasniD{T)\'ip±\'^ip± + to^r^ ^ ±_ . 

OT 

(3) 

We noted k = m/h and introduced the linear atomic den- 
sity TLiu = li/'z/p. For a dilute atomic beam, this equa- 
tion is formally equivalent to the optical equation P]). 

For an interacting atomic beam, however, the coeffi- 
cient of the nonlinear term in Eq. ^ still depends on the 
propagation parameter r. In the WKB treatment where 
the linear density reads nii3(r) = mT /p{t), this coeffi- 
cient is a constant only if no significant longitudinal po- 
tential is applied to the atomic beam, which propagates 
then with a constant longitudinal velocity. The nonlin- 
ear optical equation ([1]) and atom-optical equation ([3]) 
are then fully equivalent. 



Eq.(|4]) appears as the Schrodinger equation [i] of a 2D 
wave- function [T5| . The corresponding Hamiltonian 

H = (k + V + S , (5) 
involves three terms 

r k = f/h^ 

lv=^-fS^'Hr-r') (6) 
[ U = k^a^f^ . 

We have introduced the operator p — ~ih{Wx, Vj,) satis- 
fying [i,Pa;] = ifii associated with the transverse momen- 
tum of the atoms or photons inside the considered beam. 
The average values of these operators are defined with 
the normalized transverse wave- function ip^, following 
usual conventions in quantum mechanics. For instance, 
the second-order moment (f^) reads 

(f2) (u) = J dxdyipl {x, y, u) {x^ + y'^)ip_L {x, y, u) . 

Let us note that, because of the correspondence 2k h, 
the Hamiltonian ^ does not have the dimension of an 
energy. Accordingly, we shall refer to the three opera- 
tors K, V, U as "effective potentials" . K is proportional 
to the transverse kinetic energy, while V and U play re- 
spectively the roles of an effective 2D-contact potential 
and of an external potential. The proposed analysis of 
the propagation essentially relies on the commutation re- 
lations between these three operators. 



C. Quantum mechanical interpretation of the 
nonlinear wave equation. 

Both equations of propagation can be written in the 
generic form 

2ik^ = -eAtV'x + 7lV'±PV^± + ePa^{uy-iP^ . (4) 
ou 

with u = z, k = y^Lu/c, e = —1, 7 = Xe^^^'^/^r ^-iid 
a = P for the optical equation, and u = t, k — m/h, 

6 = 1,7 = STrniDfls and a — wj, for the atomic equation. 

One must assume that no longitudinal potential is driv- 
ing the atoms in order to obtain a nonlinear coefficient 

7 independent of the coordinate t in the atomic equa- 
tion. As stated previously, we make this assumption for 
the subsequent developments: otherwise, the equations 
of motion for the beam moments are not solved by the 
ABCD matrix approach presented here. However, the 
coefficient a of the quadratic effective potential may de- 
pend on the coordinate u in both the optical and atomic 
equations. 

If one makes the correspondence 2fc — > ft and u t, 



III. FREE NONLINEAR PROPAGATION OF 
OPTICAL AND ATOMIC BEAMS. 

We treat here the simpler case for which the parax- 
ial equation ^ contains no effective quadratic potential. 
This corresponds to the propagation of an optical beam 
in a Kerr medium of uniform linear index, or to that of an 
interacting atom laser in the absence of external poten- 
tial. We use the hidden symmetry of the paraxial equa- 
tion to derive the equations of motion for the moment 
(f^) associated with such beams. We obtain a simple 
expansion law, which we apply to discuss experimental 
results ^ on guided atom lasers. 



[i] To treat atomic beams only, one could indeed obtain directly a 
2D Schrodinger equation from Eq. JSJ instead of Eq. l|4]l, without 
having to do the correspondence 2k — > h. The great advantage 
of starting from Eq. Q is that one obtains a treatment address- 
ing simultaneously optical and atomic beams. This also shows 
the insight provided by a quantum-mechanical approach of the 
optical paraxial equation. 
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A. Effective Hamiltonian scale invariance and 
propagation of the second-order moment. 

The Hamiltonian ([5]) associated with the free propaga- 
tion equation reduces to 



HQ = e K + V . 



(7) 



The key point is to investigate how this operator trans- 
forms in a space dilatation r ^ Ar [11,] 



Ar, 7/.(r) -> A^/V(r/A), 



A2 



At + (Ar) 



To ensure the scale invariance of the Hamiltonian iJo, 
the interaction potential should transform as V^(Ar) = 
\~'^V{y). This is verified by the potential V{y) — 
l/r^, but also by the 2D-contact potential V^(r) = 
1-^(5(2) (i-) The scale invariance of Hq can be ex- 

pressed equivalently in terms of commutation with the 
generator Q of dilatations as 



[Q, Ho] = 2tHo 



(8) 



This generator reads Q — ^ (p • r + f • p), it operates at 
each point r a translation proportional to the vector r. 
The evolution of the second-order moment (identified to 
the width w'^ — (f^) for a light beam) is given by Ehren- 
fest's theorem, applied with the Hamiltonian ^ 



d(f2) 

du 



2e 



{[Q,Ho]) 



(9) 



The parabolic evolution of (f^) is indeed a direct conse- 
quence of the Hamiltonian scale invariance ([S]) : the quan- 
tity d2(f2)/du2 is proportional to the average effective 
energy (Hq) of the transverse wave- function 



2e , - > 



dw2 



which is a constant of motion. This result holds only 
for a 2D transverse space, the only dimensionality giving 
the desired scale transformation of the contact potential. 
In this case, using Eq.® and Eq.®, one obtains readily 
the free expansion law valid for optical and atomic beams 

(^')(") = ^(^o) + |(Q)o u + (f2)o . (10) 

We noted ( )o the average value of the considered op- 
erator at the initial time of the propagation, namely for 
the coordinate u = 0. 



Application to the time-of-flight analysis of an 
atom laser beam. 



the free-fall frame, or alternatively by applying Ehren- 
fest's theorem as in Eqs. ^ with an additional poten- 
tial Vg = 2k^gy (with Oy the vertical axis), one re- 
trieves a parabolic law of Eq. PH)) for the transverse width 
w2(r) — (x^) — (x)2 -I- (y2^ _ (y^2^ consider a wave- 
function corresponding to a centered atomic beam, veri- 
fying initially (f)o = and {Q)a = 0. e — I and k = m/h 
in the atomic equation, so the free atomic Hamiltonian 
is Hq ^ K + V. The width evolution reads 



w'{r)^^[{K)o + {V)o)T' + w'iO). (11) 

The quantity {K)o — / h? {/S.v'^) q reflects the contri- 
bution of the transverse kinetic energy to the beam di- 
vergence, while the quantity (V)o reflects that of mean- 
field interactions. This second contribution is usually dis- 
carded in the time-of-flight analysis of atom laser beams. 

The expansion law pTjl can be readily applied to in- 
vestigate interactions effects on the transverse dynamics 
of a guided atom laser. Assuming that one knows suf- 
ficiently the initial wave-function profile as to estimate 
(y)o, this equation gives the transverse velocity disper- 
sion (Aw2)g as a function of the width expansion w{t) 
measured experimentally 



(A«2 



2t2 V 



w'^{T)-w'^{Q)~h^/m^{V)o 



Interaction effects have been neglected in the analysis of 
the expansion given in 0, @ , in which the expansion is 
simply attributed to the initial transverse kinetic energy. 
This leads to a different estimate of the transverse veloc- 
ity dispersion {IS.v'^) free, expressed as a function of the 
beam width 



{Av'^)free 



1 

2^2 



w 



'■{0)) 



Discarding interaction effects in the beam expansion thus 
yields the following error on the transverse velocity 



(Aw2 



'free 



(Ai;2)o 



{V)o 
{K)o 



(12) 



For the quasi-monomode beam reported in 0], one can 
estimate the ratio {V)o/{K)q by modelling the initial 
transverse wave-function with a Gaussian profile verify- 
ing (r)o = (p)o = [ii]. Ea. p2)) shows then that the 
value (Aw2) ^^ee overestimates the correct square veloc- 
ity dispersion (Ai;2)q as {Av^)free = {l+2niDas) {Av'^) Q , 
thus of roughly 18% in [sj with the reported linear den- 
sity. In fact, interaction effects may indeed affect signif- 
icantly the flight of the atomic cloud, and can only be 
safely neglected in the hmit nioas <C 1 [30l |. 



The previous result can be extended easily to account 
for the presence of a linear gravitational potential in 
the propagation of an atomic beam. By considering 



[ii] To evaluate this error for a multiple-mode atomic beam Q|, one 
needs to know the mode decomposition of the transverse wave 
profile in order to estimate (V')o. 
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IV. ABCD MATRIX ANALYSIS OF THE 
NONLINEAR PROPAGATION OF OPTICAL 
AND ATOMIC BEAMS IN PRESENCE OF A 
QUADRATIC POTENTIAL. 



We now address the nonhnear propagation described 
by the fuU Eq. ^ and obtain an exact ABCD matrix 
analysis of the beam evolution. As seen previously, the 
additional effective potential contained in this equation 
reflects either a position-dependent linear permittivity of 
the form er{x, y, z) = e°(l — /3(z)^r^) for an optical beam, 
either a time-independent quadratic potential ?7^(r, z) = 
m/2 X for an atomic beam. We remind that, 

although they are assumed time-independent, the coeffi- 
cient /3 and the trapping frequency lo±_ may depend on the 
longitudinal coordinate z or parameter r. The presented 
treatment can thus address atom laser experiments fo- 
cusing the beam thanks to a time-independent quadratic 
transverse potential which sharpens along the propaga- 
tion axis. 

We proceed as follows. Using commutation relations 
between the three operators K, V, U defined in Eqs.(l6]), 
we derive the equations of motion for the second-order 
moment. This approach is inspired from the treatment 
of Ghosh [2^, but indeed these equations had been de- 
rived several times before [H, [H, [ij, HI] • We then intro- 
duce a nonlinear quality factor for matter- wave beams, 
the invariance of which is related to the hidden symme- 
try of the wave equation. This parameter enables the 
definition of a nonlinear complex radius of curvature q 
"renormalized" by the mean- field interactions. We show 
that this radius of curvature can be propagated through 
2x2 ABCD matrices identical to those used for linear 
optical and atomic beams. This conclusion agrees with 
the results obtained by Porras et al. [l^ in nonlinear 
optics. To ease the comparison between the linear and 
nonlinear regimes, an introduction to the linear ABCD 
matrix formalism for optical and atomic beams is given 
in Appendix. 

The presented treatment avoids the limitations of pre- 
vious nonlinear ABCD methods developed for matter 
waves [Slljllll- In the non-paraxial ABCD approach [sij . 
the interaction term is considered as a perturbation, as- 
sumption which becomes invalid for dense clouds or for 
long propagation times. The paraxial ABCD method 
used so far for interacting atomic beams [3l|, [111, in- 
spired from Ref. [s^l, requires the absence of external 
potential and, most critically, assumes a Gaussian-shape 
approximation of the beam likely to break down after a 
finite propagation time. The advantage of the proposed 
treatment is that it addresses interacting atomic beams 
in a fully non-perturbative manner. As in previous non- 
linear paraxial ABCD approaches [3l|, [11], this method 
applies only to guided atomic beams experiencing no sig- 
nificant longitudinal potential and propagating at a con- 
stant longitudinal velocity. This hypothesis is realistic 
for the guided atom lasers reported in the Refs. [1, [H] ■ 



A. Solutions for the propagation in a transverse 
quadratic potentiaL 

In the presence of a quadratic potential U, the Hamil- 
tonian ^ H = eK -\- V + dJ \sx\o longer scale-invariant, 
but the operators — eK + V,U,Q verify nonetheless 
the special set of commutation relations [ll| 

[Q, Ho] = 2iHo [Q, elf] = -2ieU [eU, Hq] = iia^k'^Q . 

(13) 

Thanks to this set of equations, signature of a hidden 
symmetry, the Hamiltonian can be embedded into an 
50(2,1) algebra [ll|- These relations yield, through 
Ehrenfests's theorem applied with H = Ho + ell, a closed 
set of coupled equations for the derivatives of the second- 
order moment 



d(f2) 

du 



2e 



(Q) 



= I ( (^0) - e{U) 



(14) 



dH 



du^ 



= -4a' 



:d(r 



du 



, da^ 
du 



Eqs.(fT4|) show that, for a constant coefficient a(u) = ao 
(or uniform trap frequency w_l(z) = wq), the second- 
order moment oscillates with a frequency independent 
from the strength of the nonlinearity. This property of 
nonlinear optics, established in [13], is formally identical 
to the universal low-energy mode found for 2D conden- 
sates [m . Exact solutions may be found with a periodic 
quadratic coefficient a{u) (or periodically modulated fre- 
quency uj±{t) for atomic beams) by noticing that the 
width obeys a Hill's equation (26l. [34]. 



B. Nonlinear Quality Factor and Radius of 
curvature. 

The proposed ABCD law applies to the propagation 
a parameter q which can be interpreted as a complex 
radius of curvature generalized for the nonlinear propa- 
gation. This parameter has been introduced in previous 
developments of non- linear optics [13, [3] , and here we 
extend its definition to atom optics. This definition in- 
volves an invariant of propagation, obtained from the set 
of commutation relations (jl3p behind the hidden symme- 
try, which reads 



Mf 



(15) 



For the optical propagation in which e = — 1, one re- 
trieves the nonlinear-optics quality factor [3, [S] ■ The 
connection of this optics invariant with the hidden sym- 
metry of the paraxial equation has never been pointed out 
to my knowledge. For the atomic propagation in which 
e = 1, the parameter Mj constitutes the generalization 
to paraxial interactin g be ams of the quality factor intro- 
duced by Riou et al. [221 ] for dilute atomic beams. The 
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term (_ffo) includes in M/ the additional divergence re- 
sulting from interaction effects, which were not accounted 
for in the quality factor defined in Ref. [l^l ■ The interact- 
ing quality factor of a cylindrical and centered ((r) — 0) 
atomic beam in the fundamental Gaussian mode reads 
Mf = 1 + 2niDas- In the limit of dilute atomic waves, 
one retrieves Mj = 1, so the interacting quality factor 
coincides with the parameter defined in 22]. It is also 
useful to introduce the radius of curvature R [l^, 



1 

R 



1 dw2 

2w^ du 



(16) 



With these parameters at hand, the generalized complex 
radius of curvature reads 



1 

R 



iMf 



(17) 



if and only if the ABCD matrix satisfies the differential 
equation 



^( A B 



1 

-a^{u) 



A B 
C D 



(22) 



as well as the condition {AD — BC){u) — 1 for any pa- 
rameter u. Setting the initial condition [A^ B, C, D] (0) = 
[1,0,0,1] guarantees then that Fi{u) and F2{u) coincide 
with and uP' / R at all times, thereby ensuring the 
ABCD law (UHl). The nonhnear ABCD matrix obeys 
Eg. ip^ . which is also the differential equation satisfied 
by the ABCD matrix describing the linear propagation 
of dilute atomic beams [See Appendix Eq. (jAlip ], This 
shows that in spite of the mean- field interactions, the 
usual interaction-free paraxial ABCD matrices [21] ap- 
ply rigorously to propagate the parameter q. In particu- 
lar, one retrieves the matrix 



C. Demonstration of the ABCD law for the 
nonhnear propagation. 

The previous assumptions regarding the uniform linear 
atomic density and cylindrical symmetry of the effective 
transverse potential still apply. Unless specified other- 
wise, we also assume that the beam is centered ((r) — 0) 
and that the potential contains no linear term [iii] . 

Considering the radii of curvature qi = q{0) and q2 ~ 
q{u) at different stages of the propagation, we seek to 
establish an ABCD law of the form 



92 = 



Aqi+B 
Cqi+D 



(18) 



It is sufficient that the two following relations be satis- 
fied [13 



^2 
i?2 



B ' 


2 


MfB^ 








B ' 











MfBD 



(19) 
(20) 



The two functions Fi{u) and F2(u), defined respec- 
tively by the RHS of Eqs. HH) and ([SQ]), obey the same 
first-order differential system as the functions and 
/R -obtained by considering the set of Eqs. p4)) and 
Eqs. (fT5TT6)l - 



dFi 
du 



= 2Fo 



Mf 



'^F,^'-if+F^-a^iu)F^. (21) 



A B 
C D 



1 u 
1 



for the free propagation [iv] and the matrix 



A B 
C D 



cos(q!u) a ^ sin(Q;u) 
— Q;sin(Q!u) cos(au) 



for the propagation in a quadratic index or cylindrical po- 
tential described by Eq.Q. These matrices are identical 
to those given in the Eqs. (|A2IA3[) for the linear optical 
propagation, or in the Eqs. (|A12IA13[) for the linear atom- 
optical propagation. It is a rather remarkable fact that 
the linear ABCD laws persists for an interacting beam: 
the inclusion of the invariant Mj in the definition of the 
radius of curvature q captures the nonlinear effects on the 
propagation. This is indeed a consequence of the special 
set of commutation relations resulting from the hidden 
symmetry. 

Setting u = T and a — in these matrices, one 
obtains the ABCD law for the propagation of an in- 
teracting atomic beam. An analogous property for the 
propagation of the parameter q had been established in 
a more restrictive framework: it applied only to free- 
propagating paraxial Gaussian atomic beams treated 
with a Gaussian-shape approximation (3ll . Is^ , which im- 
plies a perturbative treatment of the interaction-term. 
The result obtained here is much more general: it is non- 
perturbative, it applies to atomic beams propagating in 
transverse potentials, and it does not require any approx- 
imation once the assumptions of paraxiality, uniform lin- 
ear density of the atomic beam (and thus absence of a 



[iii] In presence of a quadratic potential, a linear gravitational po- 
tential can nonetheless be easily accounted for with a change of 
coordinates x' = x, y' = y + g/uj'^ (sag), and by considering the 
evolution of the moment {f''^) instead of {r'^). 



[iv] For the free propagation, this result can be obtained directly 
by noting that Eqs.((9ll imply (1/g)' = — and thus simply 
^(t) = q{0) + T. By replacing (r^) with in the definition JTSj 
of the quality factor, this relation can be extended to the prop- 
agation of a non-centered beam ((f) 7^ 0) in a linear potential. 
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longitudinal potential), and cylindrical symmetry of the 
transverse potential are satisfied. 

In a similar manner as in [l^, [31 , this approach could 
also be used to discuss the self-trapping in future exper- 
imental beams involving atoms with attractive interac- 
tions, such as fermions with Feshbach-tuned interactions. 
Self-trapping is expected at a critical density yielding 
Mi = 0, i.e. nuj = l/(2|as|) for Gaussian beams, be- 
yond which the beam collapses. 



APPENDIX A: THE ABCD MATRIX 
FORMALISM FOR LINEAR PARAXIAL OPTICS 
AND ATOM OPTICS. 

This Appendix gives an introduction to the ABCD 
matrix formalism for optical and atomic beams experi- 
encing a linear propagation. 

1. The ABCD formalism in linear paraxial optics. 



V. CONCLUSION. 

To summarize, I have exploited the symmetries of 
the Hamiltonian associated with the nonlinear paraxial 
wave equation - scale invariance in the free propagation, 
specific set of commutation relations in presence of a 
quadratic external potential - to extract simple propa- 
gation laws for the width of a paraxial atomic beam of 
uniform density. These results can be applied to anal- 
yse the divergence of atom laser beams in recent exper- 
iments 0, [E] and evaluate their number of transverse 
modes. A quality factor [Eq. has been defined, 

valid for paraxial interacting atomic beams propagating 
in cylindrical potentials, thereby generalizing the param- 
eter introduced in [22( | for dilute paraxial atom lasers. 
This parameter, together with a generalized radius of 
curvature, allows one to describe with linear ABCD ma- 
trices the propagation of a paraxial interacting atomic 
beam in constant, cylindrical and quadratic potentials. 
An interesting question is the robustness of the presented 
approach to a relaxation of the hypothesis of paraxiality 
and uniform density of the atomic beam. Recent the- 
oretical work [sTi on a matter- wave resonator (36| . in- 
volving a free-falling atomic cloud bouncing on curved 
atomic mirrors, shows that the width oscillations of uni- 
versal frequency persist with a good approximation in 
spite of a time-dependent linear density. This suggests 
that this treatment could be fruitfully applied to guided 
atom lasers presenting a slow variation of this parameter. 
Last, let us point out a strong connection between the 
physics of guided atom lasers and of 2D condensates, in- 
volving a similar equation. We have seen that theoretical 
tools for the latter can also be relevant for nonlinear atom 
optics. This also suggests that several effects specific of 
2D condensates, such as the breathing modes, might be 
reproduced within the transverse profile of guided atom 
laser beams. 

The author thanks David Guery-Odelin and Yann Le 
Coq for fruitful discussions, William Guerin for com- 
ments, and Steve Walborn for manuscript reading. He 
acknowledges a very rich collaboration on atom optics 
with Christian Borde, and thanks Luiz Davidovich and 
Nicim Zagury for hospitality. This work was supported 
by DGA(Contract No 0860003), by the Ecole Polytech- 
nique and by the French Ministry of Foreign Affairs 
(Lavoisier-Brsil Grant). 



The ABCD matrix formalism, discussed in various 
textbooks [37|, was originally introduced to analyse the 
light wave propagation in linear optical systems involv- 
ing a constant or piece-wise quadratic index of refrac- 
tion [l9| . We present briefly this approach for ray and 
Gaussian optics. 

One associates to a light ray a two-component vec- 
tor map involving the ray distance to the axis r(z) and 
the slope r'(z) — dr{z)/dz. In a linear optical system, 
the propagation of this vector is given by a linear input- 
output relation of the form 



r{z2) 



A 
C 



r(zi) 
r'(zi) 



(Al) 



called ABCD law, in which the matrix coefficients de- 
pend on the propagation distance Az = Z2 — zi, and sat- 
isfy the relation AD — BC =1. In a medium of uniform 
index, the ABCD matrix is simply 



(A2) 



In a graded index medium with n{r,z) = no(l — 
l/2/3'^(z)r^), the ray propagation equation yields the ma- 
trix 
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A B 
C D 



cos(/3Az) /3-isin(/3Az) 
-/3sin(/3Az) cos(/3Az) 



(A3) 



The ABCD formalism is especially well-suited to prop- 
agate Gaussian waves. Let us consider an input beam at 
the plane z = zi with a transverse profile 



i)±{x,y,zi) = exp ( 



2i 



. , 1 1 

with — = — J . 

qi Ri kwf 

(A4) 

i?i is the radius of curvature and wi the beam width 
in the direction O^, qi is called the complex radius of 
curvature. The quantity k = nQUj/c is the longitudinal 
wave-vector of the optical wave. The beam is assumed 
to propagate between the planes z = zi and z = zi = 
z\ -f Az in a linear optical system characterized by an 
ABCD matrix. The beam propagator in this system is 
directly related to the ABCD matrix elements 37i] 



K{x, X , Az) 



ik 
2^ 



exp 



2B^ 



2x'x 



Dx^) 
(A5) 
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The output beam at z = Z2 is given by the propagation 
integral 



/ + 00 
-oo 

and can be expressed as 



dx'K{x', X, Z2 - zi)ip±{x', y, zi) , 

(A6) 



kx^ 

, exp I —i 

A + B/qi ^ V 2(72 



with the complex radius of curvature (72 given by the 
ABCD law, rewritten here for convenience 
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Aqi+B 
Cqi+D 



(A7) 



This treatment can be extended to the full Hermite- 
Gauss basis 13711. 



2. The ABCD formalism in linear paraxial atom 
optics. 

Let us now apply the ABCD matrix technique to prop- 
agate dilute paraxial atomic waves [1, [2lj. The corre- 
sponding propagation equation has the form of Eq.([2l) 
without the interaction term (jiio — 0) 



ih 



At + V^{x,t) 

Zm 



(A8) 



The potential is assumed to be quadratic of the form 
Vx — muj'^i The propagator K is give n by the Van 



Vleck formula or by the ABCD formalism [15| 



K{x, x' , t) 



2h7rB 



exp 



/ im 



\2hB 



{Ax'^ ~ 2x'x + Dx^) 



(A9) 

If one sets k — m/h as previously for atomic beams, this 
expression becomes formally identical to Eq. (|A5p up to a 
sign change in the argument of the exponential. The cor- 
responding ABCD matrix describes the one-dimensional 
motion of a classical particle moving in the potential Vx ■ 
Precisely, let us consider a classical particle initially at 
the position x{0) and with a speed ^^^(O), moving in the 
potential Vx- At the instant t = t, this particle has a fi- 
nal position x{t) and a speed Vx(t). Since the equations 
of motion are linear, the final values are related to the 
initial ones by a linear map analogous to Eq. (|A7|) 



x{t) 

Vx{t) 



Air) 
C{r) 



B{t) 
D{t) 



x{Q) 

Vx{Q) 



(AlO) 



The ABCD matrix involved in this relation is precisely 
the matrix used in the propagator of Ea. (|A9p . Let us 
verify this. One can check that the distribution defined 
by Eq. (jA9p satisfies the Schrodinger equation if and only 



if the coefficients A,B,C,D satisfy the first-order differ- 
ential equation 



d 

d7 



A B 
C D 



I 

-uj^, 



A B 
C D 



(All) 



It is easy to show that the ABCD matrix defined in 
Ea. (jA10| l obeys the same equation. Besides, when the 
time T tends towards 0, the ABCD matrix defined by 
Eq. (|A10[) tends towards identity. In this limit, the dis- 
tribution defined with this matrix in Eg. ljAQD yields a 
Dirac distribution, as expected for a propagator between 
equal instants. Eq. (|A9|1 and the ABCD matrix of 
Ea. (jA10| l thus define a distribution K which satisfies 
the Schrodinger equation as well as the correct initial 
condition. This shows that K is the propagator of the 
Schrodinger equation. 

Using Ea. (|All[ ). it is straightforward to establish the 
expression of the ABCD matrix for the free motion 



A B 
C D 



1 T 
1 



and for the motion in a quadratic potential 



A B 
C D 



cos{ujj_t) a;^"'"sin(a;^r) 
— cij_Lsin(a;^r) cos(lj_lt) 



(A12) 



(A13) 



The expression of the propagator can be readily ap- 
plied to obtain, by analogy with Gaussian optics, an 
ABCD law for the propagation of Gaussian wave- 
functions. Let us consider the complex radius of curva- 
ture, decomposed as previously along its imaginary and 
complex part as 



1 



1 



2i 



We consider the initial transverse wave- function , ini- 
tially of the form 



■0±(a;,y,ri) 



1/4 



exp 



i x"^ 



(AM) 



Note that this wave-function is properly normalized to 
unity. Its evolution is given by the propagation integral 

'(p±{x,T2) = j Ax K{x,x ,T)lpA_{x ,y' ,Ti) , 

with T = T2 — Ti. Using the optical integral (|A6p . one 
obtains the wave-function at the instant T2 



'^±{x,y,T2) = - 



1/4 



1 



/W2 



: exp 



I X 

2q2 



where qi and 52 satisfy the ABCD law of Ea. (|A7[) . 

This formalism can be readily extended to the propaga- 
tion of wave- functions in a separable quadratic transverse 
potential V± = Vx + Vy. The corresponding propagator 
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is simply the product of the one-dimensional propaga- 
tors K{x,x' ,y,y' ,t) = K^{x,x',T)Ky{y,y\T) defined by 
Eq. (jA9ll . A 2 X 2 ABCD matrix is associated with each 
direction. The propagation of a Gaussian wave-function 



is formally identical: the parameters Qx and qy satisfy a 
relation of the form of Ea. (|A7p with the ABCD matrices 
associated with their respective directions. 



V'-L(a;,y,Ti) 



TTWlxWly 



■ exp 



2(?i 



exp 



2(?ly' 
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